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The continued fraction method (also known as Leaver’s method) is one of the most effective tech¬ 
niques used to determine the quasinormal modes of a black hole. For extremal black holes, however, 
the method does not work (since, in such a case, the event horizon is an irregular singular point of 
the associated wave equation). Fortunately, there exists a modified version of the method, devised 
by Onozawa et al. [Phys. Rev. D 53, 7033 (1996)], which works for neutral massless fields around an 
extremal Reissner-Nordstrom black hole. In this paper, we generalize the ideas of Onozawa et al. to 
charged massless perturbations around an extremal Reissner-Nordstrom black hole and to neutral 
massless perturbations around an extremal Kerr black hole. In particular, the existence of damped 
modes is analysed in detail. Similarities and differences between the results of the original continued 
fraction method for near extremal black holes and the results of the new continued fraction method 
for extremal black holes are discussed. Mode stability of extremal black holes is also investigated. 
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I. INTRODUCTION 

Perturbation theory is used everywhere in Physics, 
from the three-body problem of Celestial Mechanics to 
Feynman diagrams of Quantum Field Theory. In Gen¬ 
eral Relativity, in particular, linear perturbations of a 
background metric (or, similarly, linear perturbations of 
a spin-s field in a curved spacetime) satisfy a second order 
linear PDF. In some cases, e.g. if the spacetime is spher¬ 
ically symmetric (as in a Schwarzschild black hole) or, 
more generally, Petrov-type D (as in a Kerr black hole), 
the PDF can be transformed, by separation of variables, 
into a system of ODFs which is more amenable to math¬ 
ematical treatment. 

Historically, Regge and Wheeler, in the late 1950s, 
were the first to employ perturbation theory in Black 
Hole Physics. In their seminal paper [1], after appropri¬ 
ate boundary conditions were imposed to the perturba¬ 
tion equations, the stability of Schwarzschild black holes 
was established. Later, in 1970, Vishveshwara [5] discov¬ 
ered a special class of solutions for the perturbation equa¬ 
tions of Schwarzschild black holes which are characterized 
by a peculiar set of boundary conditions: purely outgo¬ 
ing waves at spatial infinity and purely ingoing waves at 
the event horizon. These solutions, known as quasinor¬ 
mal modes (QNMs), correspond to stable modes which 
decay in time. Assuming a time dependence of 
these modes are characterized by complex frequencies w 
(the quasinormal frequencies) whose imaginary part is 
negative. In general, only a discrete set of quasinormal 
frequencies will be allowed. 

Since the discovery of QNMs, several methods have 
been devised to determine the quasinormal frequencies of 
a black hole system. The simplest one consists in approx¬ 
imating the black hole potential by a Poschl-Teller po¬ 
tential whose analytical solutions are known . Another 
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approximation technique is the WKB method [4], which 
is equivalent to finding the poles of the transmission co¬ 
efficient of a tunnelling problem in Quantum Mechan¬ 
ics. One can also apply a shooting method to match the 
asymptotic solutions at some intermediate point, as done 
by Chandrasekhar and Detweiler [^. Another important 
method is the direct integration of the perturbation equa¬ 
tion in the time domain using light-cone coordinates [^. 
The most accurate method, however, is the continued 
fraction method developed in 1985 by Leaver and 
later improved by Nollert [3]. For a detailed review on 
the available methods, we recommend [10] (see also [TTl - 
HSj for other reviews about QNMs). 

Leaver’s method was inspired by a technique due to 
Jaffe to calculate the energy eigenvalues of the H 2 ion. 
It consists in using the Frobenius method of differential 
equations to write the solution of the perturbation equa¬ 
tion as a power series around the event horizon. It can 
be shown that this series will satisfy the QNM boundary 
conditions only if a certain equation involving an infinite 
continued fraction is also satisfied. By solving the contin¬ 
ued fraction equation with a root hnding algorithm, one 
is then able to determine the QNMs. Fven though Leaver 
originally applied his method only to massless perturba¬ 
tions of Schwarzschild and Kerr black holes, it can be 
implemented in a variety of other situations unnsj. 

In this paper, we are particularly interested in QNMs 
of extremal black holes. Unfortunately, Leaver’s original 
method does not work in such cases, and the reason lies in 
the fact that the event horizon of an extremal black hole 
is an irregular singular point of the perturbation equa¬ 
tion. Therefore, unlike the case of non-extremal black 
holes (for which the event horizon is a regular singular 
point), one does not expect a non-zero radius of con¬ 
vergence of the associated power series around the event 
horizon. In 1996, however, Onozawa et al. m success¬ 
fully devised a modihcation of Leaver’s original method 
for an extremal RN black hole which does not rely on 
expansions around the event horizon. Instead, they ex- 
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pand the solution around an ordinary point of the differ¬ 
ential equation. Following the ideas of Onozawa et al., 
we show in this paper (Secs. [lT| and III) how to mod¬ 
ify the original continued fraction method to calculate 
the QNMs of an extremal Kerr black hole. We also ap¬ 
ply the method to spin-0 and spin-1/2 charged pertur¬ 
bations around an extremal RN black hole, generalizing 
the results of Onozawa et al. for non-neutral fields. Our 
numerical results, including a detailed analysis of the so- 
called damped modes for extremal Kerr black holes and 
a mode stability investigation, are presented in Sec. m 
followed by our final remarks in Sec. |V| 


II. FIELD DYNAMICS 
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^sin cos^ 9 — 2aujs cos 9 
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where A is a separation constant, K = u}(r^ -b a^) — am, 
and X = Aisuir + 2maio — a^uj^ — A. 

The description of a charged black hole system is sim¬ 
ilar. The RN line element is given by 


According to the black hole uniqueness theorems [m 
[IS, under reasonable conditions, the Kerr metric is the 
only stationary axisymmetric black hole solution of the 
vacuum Einstein field equations while the RN metric is 
the only spherically symmetric black hole solution of the 
Einstein-Maxwell field equations. The first one describes 
the spacetime of a rotating black hole with mass M and 
specific angular momentum a < M, and the second one 
describes an electrically charged black hole of mass M 
and charge Q < M (throughout this paper we use units 
in which G = c = I). 

The Kerr line element, in Boyer-Lindquist coordinates 
{t,r,9,(j)), is given by 


ds^ = - 1- 


2Mr\ 

) 


, 2 AMrasw? 0 ,, , , 2 

dt^ - 5 - dcbdt + -^dr^ 

A 


+p^d9^ 


-b 


2Mra^ sin^ 9 


sin^ 9 d(j)^, 

( 2 . 1 ) 


where p^ = -b a? cos^ 9 and A = r^ — 2Mr -b a^. 
If a < M, the roots of the A function, namely = 
M + — 0 ? and r- = M — y/M'^ a^, correspond 

respectively to the event horizon and the Cauchy surface 
of the black hole. If a = M, these two surfaces merge 
into the surface r_|_ = r_ = M, which corresponds to 
the event horizon of an extremal Kerr black hole. Fi¬ 
nally, if a < M, there is no event horizon and the metric 
corresponds to a naked singularity at r = 0. 

In the early 1970s, Teukolsky derived one of the most 
notable facts about the Kerr metric: that its pertur¬ 
bations are separable (in the frequency domain) when 
the Newman-Penrose formalism is employed |17j . In 
fact, let '0 describe any neutral, massless spin-s field 
(|s| = 0,1/2,1,2) and write the corresponding wave 
equation in the Kerr background using the ansatz ip = 
R(r)S'(0)e™'^““C As a result, the wave equation re¬ 
duces to the so-called Teukolsky equations [171 [TB] for 


A 


= - ‘^dt^ + ^dr^ + {dd"^ + 


sm 




(2.4) 


where A = r'^ — 2Mr + Q^. Its event horizon is located at 
r+ = M -b y/which reduces to r+ = M in the 
extremal case {Q = M). Furthermore, the study of spin-0 
and spin-1/2 electrically charged, massless fields around a 
RN black hole can be simplified if the same ansatz used 
for Kerr perturbations is employed. Indeed, by doing 
so, one can show that the Klein-Gordon and the Dirac 
equations, in the RN metric, can be separated and recast 
as new equations which are analogous to the Teukolsky 
equat ions [T^ 12(1] . The angular equation assumes the 
form \2.2\ with aw = 0, while the radial one assumes the 
form ( |2.3[ ) with K = ujr'^—qQr and X = 4isujr—2isqQ—X 
{q denotes the charge of the held). 


In order to solve the perturbation equations (2.2) and 


(2.3), one needs appropriate boundary conditions. For 


the angular function, it is natural to require regularity 
at the singular points 9 = 0 and 9 = n, transforming 
eq. (2.2) into a Sturm-Liouville problem for the separa¬ 


tion constant. For a given set of parameters aio, s and 
m, only a discrete set (indexed by the new parameter £) 
of constants A is allowed. The associated functions are 
called the spin-weighted spheroidal harmonics [HI [2T] . 
In the particular case of aoj = 0, valid for the RN met¬ 
ric, these angular functions reduce to the spin-weighted 
spherical harmonics, for which the separation constant 
can be determined analytically as A = (£ — s){i -b s -b I). 
Additionally, the parameters satisfy the following con¬ 
straints: m is an integer (bosonic case) or a half-integer 
(fermionic case), (. is & non-negative integer (bosonic 
case) or a non-negative half-integer (fermionic case), 
I > |s|, and —l<m<l. 

Regarding the radial equation, on the other hand, we 
impose the usual boundary conditions for QNMs: purely 
ingoing modes at the event horizon (since no classical per¬ 
turbation can escape the black hole) and purely outgoing 
modes in the asymptotic limit r —> 00 . In these limits, 
eq. (2.3) can be solved analytically, yielding the following 
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boundary conditions (for an extremal black hole): 

R{r) (X {r — MY^ , when r —>■ M, (2.5) 

Rir) (X when r ^ oo, (2.6) 

where Jq = iM{2Moj — m), Ji = —2s — 2iMuj, and 
J 2 = — 1 + AiMoj in the case of an extremal Kerr black 
hole, and Jg = iM{Muj — qQ), Ji = —2s — 2iMuj + iqQ, 
and J 2 = — 1 + AiMuj — 2iqQ for an extremal RN black 
hole. 


III. CONTINUED FRACTION METHOD 


As explained in the introduction, the original contin¬ 
ued fraction method devised by Leaver, which consists in 
expressing the solution of the radial Teukolsky equation 
as a power series around the event horizon, works only 
for non-extremal black holes. Therefore, inspired by the 
ideas of Ref. m for neutral perturbations around an ex¬ 
tremal RN black hole, we expand the solutions of the 


radial equation (2.3) around the ordinary point r = 2M. 


More precisely, we write 


R{r) = e^‘^^e^{r-MY^r^^^an( / , 

n^O ^ 

(3.1) 

where the expansion coefficients a„ satisfy a five-term 
recurrence relation (the last equation holds for n > 3): 

0102 -I- /3iai -I- 7 iao = 0, (3.2) 

0:203 -I- ,0202 -I- 7201 -I- (5200 = 0, (3.3) 

OnOn+l T /5nOn -f 7nOn—1 T ^n^n—2 T ^n^n—3 0- 

(3.4) 

The coefficients an, Pn, In, <5™ and e„, for both Kerr 
and RN extremal black holes, are given explicitly in the 
Appendix. 


Note that the boundary conditions (2.5) and (2.6) are 


automatically satisfied provided that the power series is 
everywhere convergent. One possible way to test con- 
vergen ce is to analyse the large-n behaviour of a„_|_i/a„. 
Using (3.4), it is possible to find four different asymptotic 


solutions, namely 


and 


Ort+l 


On-t-1 


= 1±2 


= - 1±2 


—iMoj 


n 


+ 0{n-^), 


(3.5) 


+ (“) 


Since |a„+i/a„| —>■ 1 in all fou r cases, the ratio test guar¬ 
antees that the sum in (3.1) is convergent as long as 
M < r < oo. However, convergence at r = M and at 
r = oo, i.e. absolute convergence of the sum is also 


needed to ensure the fulfilment of the QNM boundary 
conditions. With the help of Raabe’s test, we find that 


the solutions associated with the minus sign in (3.5) or 


the plus sign in (3.6) are compatible with QNMs. The 


only exception happens when iMuj or Jq are positive real 
numbers, since then the coefficient of the term be¬ 

comes purely imaginary [^. 

Equation (3.4), being a five-term linear recurrence rela¬ 


tion, possesses four independent solutions, each one yield¬ 
ing one of the four possible asymptotic behaviours dis¬ 
cussed above. (In other words, one has the freedom to 
choose Oo, ai, 02 and 03 , from which all the other a^s will 
be determined.) Futhermore, equations (3.2) and (3.3) 
can be thought as boundary conditions for (3.4), reduc¬ 


ing the degrees of freedom from 4 to 2. Equivalently, one 
can perform two successive Gaussian eliminations in (3.4) 


to transform eqs. (3.2)-(3.4) into a three-term recurrence 
relation for a„: 


a"a„+i -I- /3"a„ -|- 7"o„_i =0, n > 1, 


(3.7) 


where the new recurrence coefficients a", /3" and 7 " 
are obtained recursively from the original ones [see 
eqs. (A.20) and (A.21) in the Appendix]. 

For a given frequency uj, this three-term recurrence 
relation possesses two linearly independent solutions. 
In general, these solutions will produce divergent sums 
o-n and, therefore, will be incompatible with the QNM 
boundary conditions. Consequently, in order to deter¬ 
mine the QNMs, we need to find the frequencies uj for 
which one of the solutions of ( |3.7[ ) will satisfy the con¬ 
vergence condition discussed above (and related to the 
asymptotic behaviour of a„+i/a„). 

Note that the absolute convergence of the series ^ a„ 
is equivalent to the convergence of the sums 'YjC-n '■= 
J2o-2n and J2dn ■= ^a 2 n+i- By rearranging and ma¬ 
nipulating the three-term recurrence relation (3.7) for 
three successive integers, one can eliminate all even- 
numbered (or, similarly, all odd-numbered) terms from 
the equations to decouple c„ from ci„, obtaining 

UnCn+l + PlCn + 7 nCn-l =0, n > 1, (3.8) 

andn+i+Pldn + lldn-i = Q, n>l, (3.9) 

where the new recurrence coefficients are given in terms 
of the original ones in eqs. ( A.25[ ) and (A.27). See the 
Appendix for all the details of the decoupling process. 
We remark that the decoupling in Ref. El assumes that 
the coefficient 5n is zero, while our decoupling works for 
(5„ 7 ^ 0. 

We now invoke the theory of three-term recurrence re¬ 
lations, and its relation to continued fractions [52], to de¬ 
termine the quasinormal frequencies. More precisely, we 
use the fact that On is a minimal solution of (3.7) if, and 
only if, the infinite continued fraction ■ ■ ■ 

Pi P2 P3 

converges. Furthermore, in case of convergence, the min¬ 
imal solution satisfies 

(H ^ 7i a'iija'iii 
ao P'{- P'^- P'i- ■ ■ ■ 


(3.10) 
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Of course, similar expressions will hold for Ci/cq and 
di/do provided that c„ and dn are, res pecti vely, mini - 
mal solutions of the recurrence relations ( |3.8[ ) and (3.91. 
As noted by Leaver for non-extremal black holes, mini¬ 
mal solutions are exactly the ones which correspond to 
the QNM boundary conditions. 

These continued fractions are not all independent; they 
can be re lated through (3.21. Mor e prec isely, after divid¬ 
ing (3.2) by Co, we use relation (3.10) above (and the 


analogous expression for Ci/cq) to obtain 


By invoking the fact that the sum in (3.13) is convergent 
only if bn is a minimal solution of (3.15), we find 


bi 

bo 


/?o 


Qn 


_7i_«i72 «|73 

/??- 0^2- Pl- 


(3.19) 


which, together with (3.11), can be used to determine the 


QNMs of an extremal Kerr black hole. 


IV. NUMERICAL RESULTS 


ai -hpi-l-7i = -Qfr 

Co Oo 


-Pi 


7i a'iija'iii 

P'i- P'i- P'i- 


Tl 72^273 

P°2- PI- ■ ■ ■ 

+ 71 = 0. 


(3.11) 


Alternatively, we could have used (3.2) and (3.3) to ob¬ 
tain 


72 _ ^2+ P2^ 

Pi do Pi 'yi+ai'^ 


= 0 , 


(3.12) 


which relates the even and odd numbered sequences after 
we substitute ci/co and di/do by their associated con¬ 
tinued fractions. Both equations (3.11) and (3.12) are 


equivalent and, therefore, the results of any of them can 
be used as a consistency check for the results obtained 
with the other one. 

For an extremal RN black hole, since the separation 
constant A can be determined analytically, one can di¬ 
rectly solve the continued fraction equation (3.11) to de¬ 
termine the quasinormal modes. However, for an ex¬ 
tremal Kerr black hole, since not only w, but also the 
separation constant A is not known a priori, we also need 
to solve the angular Teukolsky equation numerically. The 
mathematical treatment is exactly the same one used for 
non-extremal black holes [7] : we define u = cos 9 and ex¬ 
pand the solution Siu) around the regular singular point 
u = —1 using the ansatz 


S'(m) = e“‘^“(l - ii)T ^ (1 + u) 


\n+ 


|m-s| 


(3.13) 


n—0 


where the expansion coefficients satisfy the following 
three-term recurrence relation: 


We now proceed to employ the technique described 
in the previous section to determine the QNMs of ex¬ 
tremal black holes. For a RN black hole, we directly solve 
eq. (3.11) for the unknown ut. For a Kerr black hole, on 


the other hand, equations (3.11) and (3.19) form a system 
of coupled, algebraic equations for A and lo and, there¬ 
fore, need to be solved simultaneously. In both cases, 
we truncate the continued fractions at some sufficiently 
large order and use a root-finding algorithm to solve the 
corresponding equation(s) for the unknown(s). We then 
increase the number of terms in the continued fractions 
and repeat the procedure until the desired accuracy is 
attained. We also compare our results for extremal black 
holes with the quasinormal frequencies of near extremal 
black holes {a/M = 0.999 or Q/M — 0.999), which we 
calculate using Leaver’s original method. Our imple¬ 
mentation of Leaver’s method follows the steps discussed 
in [7] for a Kerr black hole and in [23] for a RN black hole, 
and, therefore, is not repeated here. In order to test it, 
we have used the tabulated values of Refs. jT] |53J 1^ . 
finding excellent agreement between them and our im¬ 
plementation. 

As explained in Refs. [3 uni mi, the spheroidal har¬ 
monics are invariant under the transformation (s —> —s, 
A —> A -I- 2s). This means that the Kerr black hole 
system is symmetric with respect to the transformation 
(to ——TO, u) —>■ —w*, A —>■ A*). The RN black hole 
system, on the other hand, is symmetric with respect to 
the transformation {q —>■ —q, uj —to*). Consequently, 
without loss of generality, we can assume (unless other¬ 
wise stated) that s is non-positive and that Re{Muj) > 0. 


A. Extremal Kerr black holes 


ccq^i + Pobo — 0, (3-14) 

anbn+i + Pibn + = 0, n > 1, (3.15) 

with Of®, Pn and 7 ® given by 

a® = —4(1-I-n)(l-I-n-I- |to — s|); (3.16) 

Pn = -I- |to — s| (—4aa; -I- |to -I- s| -I- 2n -I- 1) 

— 2{2aoj -I- A) — Aaoj{2n -I- s) -I- (2n -|- 1 ) |to -|- s| 

-I- mp 2n{n -|- 1) — s{s + 2); (3-17) 

7 ® = 2auj [2(n -|- s) -|- |to — s| -I- |to -I- s|]. (3.18) 


We start our analysis with to = 0 perturbations around 
a Kerr black hole. The quasinormal frequencies of an ex¬ 
tremal black hole, together with the quasinormal frequen¬ 
cies of a near extremal black hole, are presented in table 
|l| The associated separation constants A are also shown. 
Note that the relative difference between the results is 
minimal: for Re{Muj), it ranges from 0.008% (s = —1, 
£ = 2, n = 1) to 0.04% (s = £ — 0, n = 1), and for 
Im{Muj) it ranges from 0.004% (s = £ = 0, n = 1) to 
0.13% (s = —2, £ = 2, n = 0). On the one hand, our 
results show that the use of Leaver’s method for near ex¬ 
tremal Kerr black holes is an excellent way to estimate 
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TABLE I. Comparison of the QNM frequencies Muj and the associated separation constants A between a near extremal 
{a/M = 0.999) and an extremal Kerr black hole for m = 0. For each set of parameters s and i, we have calculated the least 
damped mode (n = 0) and the first overtone (n = 1). 


s 

£ 

n 


a/M = 

0.999 



a/M 

= 1 





Re(Ma;) 

lm{Muj) 

Re(A) 

Im(A) 

Re{Muj) 

lm{Muj) 

Re(A) 

Im(A) 

0 

0 

0 

0.110265 

-0.089439 

-0.001378 0.006564 

0.110245 

-0.089433 

-0.001380 0.006575 

0 

0 

1 

0.062498 

-0.318852 

0.032404 

0.013144 

0.062473 

-0.318840 

0.032468 

0.013164 

0 

1 

0 

0.314946 

-0.081771 

1.944568 

0.030907 

0.314986 

-0.081714 

1.944436 

0.030952 

0 

1 

1 

0.281435 

-0.253809 

1.991284 

0.085574 

0.281392 

-0.253686 

1.991243 

0.085691 

-1 

1 

0 

0.274777 

-0.075305 

1.972061 

0.016634 

0.274828 

-0.075232 

1.971989 

0.016655 

-1 

1 

1 

0.240106 

-0.234599 

1.999206 

0.044984 

0.240053 

-0.234445 

1.999185 

0.045035 

-1 

2 

0 

0.500902 

-0.079434 

5.907613 

0.029813 

0.501013 

-0.079365 

5.907384 

0.029852 

-1 

2 

1 

0.479397 

-0.241175 

5.934462 

0.087017 

0.479437 

-0.240991 

5.934283 

0.087130 

-2 

2 

0 

0.424998 

-0.071899 

3.907902 

0.032246 

0.425145 

-0.071806 

3.907644 

0.032281 

-2 

2 

1 

0.402689 

-0.218528 

3.940442 

0.092568 

0.402744 

-0.218283 

3.940242 

0.092663 

-2 

3 

0 

0.664945 

-0.076814 

9.856062 

0.033400 

0.665132 

-0.076735 

9.855691 

0.033441 

-2 

3 

1 

0.649315 

-0.231889 

9.877945 

0.098717 

0.649443 

-0.231668 

9.877614 

0.098836 


the actual values for the extremal case. On the other 
hand, our method provides the most accurate technique 
available for the determination of QNMs of extremal Kerr 
black holes. 

We now turn our attention to Kerr QNMs with m > 0 
(corotating modes). As first predicted by Detweiler 
in [25] after an analytical treatment of the near extremal 
regime, when the black hole approaches extremality, its 
quasinormal frequencies approach the purely real value 
Muj = ml2. This behaviour, further studied in Refs. pSk 
E], suggests the existence of QNMs with exactly van¬ 
ishing imaginary parts for extremal Kerr black holes. 
Ideally, we would like to apply the continued fraction 
method to determine if these QNMs modes with van¬ 
ishing imaginary parts are indeed part of the quasinor¬ 
mal spectrum of extremal black holes. However, when 
Muj = ml2, there is no guarantee that the meth od 
will work, since we have Jq = 0 and, therefore, (2.5) is 
not correct anymore. Nevertheless, as we increase the 
number of terms in the continued fractions, the solu¬ 
tions of eqs. (3.11) and (3.19) seem to indeed approach 
Muj = m/2. 

Besides these zero-damping modes, near extremal Kerr 
black holes exhibit, for some set of parameters i and m, 
QNMs whose imaginary part approaches a non-zero value 
as ajM —> 1 [5DH5^ . These modes, referred to as damped 
modes, can be obtained for extremal black holes using 
the method discussed in this paper. In particular, the 
damped modes we obtain for s = —2, £ = 2, m = 1 and 
s = 0, £ = 10, m = 7, presented in Fig. are compatible 
with those obtained in Figs. 8 and 9 of Ref. [31] for near 
extremal black holes. 

We also analyse the damped modes for s = 0, £ = 2, 
m = 1, which were discussed in Appendix D of Ref. [31j 
due to their peculiar behaviour. Our results, shown in 
Fig. [^together with the zero-damped modes, are in agree¬ 
ment with Fig. 14 of Ref. |31j . In particular, we were able 
to resolve the spiralling behaviour of the third damped 
mode as a/M is increased. 

The condition for the existence of such damped modes 


is given by 0 < m/(£ -I- 1/2) < 0.744 and, even though 
this criteria was found using the eikonal limit £ ^ | s |, it 
has been shown to be accurate even for low £. Recently, 
however, Hod [33] argued that these modes should be 
present even when this condition is not satisfied. Zim¬ 
merman et al. |34| . using numerical techniques, were not 
able to find any of these modes when the condition does 
not hold. 


In our numerical simulations, we chose to study not 
only the scalar and gravitational perturbations consid¬ 
ered in Refs. (3011311, but also Dirac and electromagnetic 
fields. By varying £ and m, we search for damped modes 
of extremal Kerr black holes and perform, for the first 
time, a detailed analysis of their behaviour as the pa¬ 
rameters change. Our results are presented in tables [H]- 
[v] for those pairs (£, m) for which damped modes were 
found, the fundamental frequency is calculated. These 
frequencies, and the associated separation constants, are 
plotted in Fig. In general, for fixed m, when £ in¬ 
creases, the fundamental damped modes oscillate more 
rapidly and become less stable. The real and imaginary 
parts of the associated separation constants A increase 
[following closely the spin-weighted spherical harmonics 
eigenvalues (£—s)(£-|-s-|-l)]. On the other hand, for fixed 
£, when m increases, the fundamental damped modes also 
oscillate more rapidly, but become more stable. The real 
and imaginary parts of the separation constants A both 
decrease. 


Our findings are summarized in Fig. [^ where we indi¬ 
cate which triplets (s, £, m) allow the existence of damped 
modes. In particular, for s = 0 and s = —2, our re¬ 
sults are in complete agreement with Fig. I of Ref. m- 
Therefore, regarding the controversy surrounding the un¬ 
restricted existence of damped modes, our results are 
in agreement with Refs. (30] [33]. In particular, fol¬ 
lowing the discussion in Ref. (34], we have searched 
around Muj = to/ 2 -£ (0.162 - i0. 035)e~ ^-^^^" for pos¬ 
sible solutions of eqs. (3.11) and (3.19) which violate 
0 < to/(£ -I- 1/2) < 0.744 when s = —2, £ = m = 2, 
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Re(Mw) 

FIG. 1. (Colors online.) Parametric plots of the real and 
imaginary parts of the first damped modes for (s = —2, i = 2, 
m = 1) (upper plots) and (s = 0, £ — 10, m = 7) (lower 
plots). All lines start at a/M = 0.9, end at a/M = 0.99999, 
and are obtained using Leaver’s original method; red marks, 
on the other hand, correspond to extremal black holes and 
are obtained using the method discussed in this paper. 


but no damped QNM solution was found. 

Finally, for m < 0 (counter-rotating modes), only 
damped modes are allowed. Once again, we compare 
our numerical results with the QNMs of a near extremal 
Kerr black hole for several parameters s, £, and m, as 
shown in table |VI[ Similarly to the m = 0 case, we find 
excellent agreement between the results. 


B. Extremal RN black holes 

The first (and, until now, only) work to discuss the 
implementation of a continued fraction technique to ex¬ 
tremal black holes was Ref. m. As explained before, 
they considered only neutral perturbations of extremal 
RN black holes and their implementation relied on the 
fact that the coefficient 6„ in eq. (3.4) vanishes for scalar, 



electromagnetic, and gravitational fields (spin-1/2 per¬ 
turbations were not considered). Our analysis for ex¬ 
tremal RN black holes, on the other hand, considers spin- 
0 and spin-1/2 charged perturbations. Furthermore, our 


FIG. 2. (Golors online.) Upper plot is a parametric plot of 
the real and imaginary parts of the first damped and zero¬ 
damping modes for (s = 0, ^ = 2, m = 1). All lines start 
at a/M = 0.95, end at a/M = 0.99999, and are obtained 
using Leaver’s original method. Note that, as extremality is 
approached, the zero-damping modes converge towards the 
purely real value m/2, while the damped modes approach 
complex numbers with non-zero imaginary part. The damped 
modes at extremality, indicated by red dots, are obtained with 
the method described in this paper. The bottom plot is a 
close-up view of the third damped mode (and its spiralling 
behaviour). 


technique works when Sn ^ 0, which is always the case 
for charged perturbations around an extremal black hole 
(and also for neutral perturbations around an extremal 
Kerr black hole). 

Contrary to Kerr QNMs, RN QNMs do not depend 
on the azimuthal number m. They depend, however, 
on an analogous (dimensionless) quantity, gQ, which ac¬ 
counts for the electromagnetic interaction between the 
black hole and the perturbation field. Additionally, while 
m is a discrete parameter, qQ is continuous. Flaving this 
in mind, we now analyse the behaviour of the quasinor¬ 
mal frequencies for an extremal RN black hole as the 
electromagnetic interaction term is varied. 

In our simulations, we start from \qQ\ 1 and track 
the quasinormal modes as qQ increases. Our results, pre¬ 
sented in tables IVTll and fVllll for several combinations of 
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FIG. 3. (Colors online.) Real and imaginary parts (left plots) of the frequencies of the most stable damped modes, and 
the associated separation constants (right plots), for scalar, Dirac, electromagnetic and gravitational perturbations. Note the 
regular behaviour of the QNMs as the parameters I and m are varied. 












































































TABLE II. Damped modes for scalar perturbations of an extremal Kerr black hole. For a given set of parameters i and m, the 
fnndamental frequency is given. Asterisks indicate pairs {i, m) for which no damped modes were fonnd. 


m\£ 

0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

0 

0.110245 
-0.0894331 

0.314986 
-0.0817141 

0.524122 
-0.0813231 

0.733303 

-0.0811681 

0.942564 

-0.0811041 

1.151863 

-0.0810711 

1.361184 

-0.0810521 

1.570518 
-0.0810401 

1.779860 

-0.0810321 

1.989208 

-0.0810261 

2.198560 

-0.0810221 

1 

X 

* 

* 

0.664311 
-0.0560541 

0.861758 

-0.0660051 

1.065265 

-0.0703481 

1.271142 

-0.0727611 

1.478194 

-0.0742891 

1.685913 

-0.0753431 

1.894048 

-0.0761121 

2.102459 

-0.0766981 

2.311064 

-0.0771591 

2 

X 

X 

* 

* 

1.071595 

-0.0322381 

1.239028 

-0.0506081 

1.427553 

-0.0591981 

1.624331 

-0.0640901 

1.825248 

-0.0672241 

2.028545 

-0.0693951 

2.233339 

-0.0709841 

2.439136 

-0.072195 

3 

X 

X 

X 

* 

* 

1.520325 

-0.0134961 

1.647474 

-0.0363971 

1.814968 
-0.0483861 

1.998597 
-0.0554881 

2.190465 

-0.0601301 

2.387137 

-0.0633811 

2.586865 

-0.0657771 

4 

X 

X 

X 

X 

* 

* 

2.000665 

-0.0015981 

2.081394 

-0.0239281 

2.223730 

-0.0382531 

2.391239 

-0.0471821 

2.571813 

-0.0531631 

2.760105 

-0.0574141 

5 

X 

X 

X 

X 

X 

* 

* 

* 

* 

2.536729 
-0.0135171 

2.650983 

-0.0289891 

2.800129 

-0.0393071 

2.967644 

-0.0464261 

6 

X 

X 

X 

X 

X 

X 

* 

* 

* 

* 

3.010459 
-0.0054761 

3.094502 
-0.0207131 

3.223556 

-0.0319491 

7 

X 

X 

X 

X 

X 

X 

X 

* 

* 

* 

* 

3.500258 

-0.0004701 

3.552506 

-0.0135231 

8 

X 

X 

X 

X 

X 

X 

X 

X 

* 

* 

* 

* 

* 

* 

9 

X 

X 

X 

X 

X 

X 

X 

X 

X 

* 

* 

* 

* 

10 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

* 

* 


TABLE III. Damped modes for Dirac perturbations of an extremal Kerr black hole. For a given set of parameters I and m, 
the fundamental frequency is given. Asterisks indicate pairs (£, m) for which no damped modes were fonnd. 


7n\i 

1/2 

3/2 

5/2 

7/2 

9/2 

11/2 

13/2 

15/2 

17/2 

19/2 

21/2 

1/2 

* 

* 

0.473956 
-0.0679991 

0.682532 

-0.0733841 

0.891515 

-0.0756431 

1.100653 

-0.0768731 

1.309870 

-0.0776461 

1.519131 

-0.0781751 

1.728421 

-0.0785611 

1.937729 

-0.0788541 

2.147050 

-0.0790851 

2.356380 

-0.0792701 

3/2 

X 

* 

* 

0.857821 

-0.0427221 

1.042609 
-0.0578421 

1.240161 

-0.0645631 

1.442531 

-0.0683121 

1.647263 

-0.0706881 

1.853332 

-0.0723251 

2.060233 

-0.0735171 

2.267687 

-0.0744241 

2.475525 

-0.0751371 

5/2 

X 

X 

* 

* 

1.289318 

-0.0214091 

1.437287 

-0.0429821 

1.616147 

-0.0535411 

1.807036 

-0.0596471 

2.003966 

-0.0635891 

2.204377 

-0.0663301 

2.406983 

-0.0683391 

2.611063 

-0.0698731 

7/2 

X 

X 

X 

* 

* 

1.756205 

-0.0060651 

1.859807 

-0.0296371 

2.015114 

-0.0430541 

2.191116 

-0.0511801 

2.377718 

-0.0565471 

2.570523 

-0.0603311 

2.767286 

-0.0631291 

9/2 

X 

X 

X 

X 

* 

* 

* 

* 

2.305459 

-0.0182091 

2.433825 

-0.0333531 

2.592399 

-0.0430831 

2.766705 

-0.0496901 

2.950337 

-0.0544211 

11/2 

X 

X 

X 

X 

X 

* 

* 

* 

* 

2.770804 

-0.0089871 

2.869788 

-0.0245861 

3.008991 

-0.0354651 

3.169438 

-0.0430981 

13/2 

X 

X 

X 

X 

X 

X 

* 

* 

* 

* 

3.253231 

-0.0023921 

3.321032 

-0.0168581 

3.439352 

-0.0283981 

15/2 

X 

X 

X 

X 

X 

X 

X 

* 

* 

* 

* 

* 

* 

3.785962 

-0.0102691 

17/2 

X 

X 

X 

X 

X 

X 

X 

X 

* 

* 

* 

* 

* 

* 

19/2 

X 

X 

X 

X 

X 

X 

X 

X 

X 

* 

* 

* 

* 

21/2 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

* 

* 
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TABLE IV. Damped modes for electromagnetic perturbations of an extremal Kerr black hole. For a given set of parameters I 
and m, the fundamental frequency is given. Asterisks indicate pairs {I, m) for which no damped modes were found. 


m\i 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

0 

0.274828 
-0.0752321 

0.501013 
-0.0793651 

0.716936 
-0.0801971 

0.929877 

-0.0805231 

1.141500 

-0.0806841 

1.352423 

-0.0807761 

1.562930 

-0.0808331 

1.773167 

-0.0808711 

1.983221 

-0.0808971 

2.193145 

-0.0809161 

2.402972 

-0.0809311 

1 

* 

* 

0.642174 

-0.0517541 

0.845425 
-0.0643821 

1.052494 

-0.0694971 

1.260686 

-0.0722361 

1.469351 
-0.0739341 

1.678255 

-0.0750861 

1.887296 

-0.0759181 

2.096423 

-0.0765461 

2.305606 

-0.0770371 

2.514828 

-0.0774311 

2 

X 

* 

* 

1.059445 

-0.0288711 

1.227322 

-0.0491911 

1.417480 

-0.0584271 

1.615647 

-0.0636061 

1.817659 

-0.0668921 

2.021821 

-0.0691531 

2.227310 

-0.0708001 

2.433676 

-0.0720501 

2.640651 

-0.0730301 

3 

X 

X 

* 

* 

1.514822 

-0.0110711 

1.639194 

-0.0351601 

1.807018 
-0.0476831 

1.991365 

-0.0550381 

2.183932 

-0.0598171 

2.381218 

-0.0631521 

2.581470 

-0.0656011 

2.783713 

-0.0674721 

4 

X 

X 

X 

* 

* 

* 

* 

2.075779 

-0.0228721 

2.217507 

-0.0376141 

2.385224 

-0.0467621 

2.566187 

-0.0528681 

2.754888 

-0.0571961 

2.948523 

-0.0604081 

5 

X 

X 

X 

X 

* 

* 

* 

* 

2.533259 

-0.0126591 

2.646200 

-0.0284141 

2.795154 

-0.0389171 

2.962808 

-0.0461481 

3.141538 

-0.0513641 

6 

X 

X 

X 

X 

X 

* 

* 

* 

* 

3.008760 

-0.0048541 

3.090940 

-0.0202051 

3.219481 

-0.0315891 

3.372458 

-0.0397231 

7 

X 

X 

X 

X 

X 

X 

* 

* 

* 

* 

* 

* 

3.549997 

-0.0130871 

3.656826 

-0.0248391 

8 

X 

X 

X 

X 

X 

X 

X 

* 

* 

* 

* 

*■ 

* 

4.021952 

-0.0071761 

9 

X 

X 

X 

X 

X 

X 

X 

X 

* 

* 

*■ 

* 

* 

* 

10 

X 

X 

X 

X 

X 

X 

X 

X 

X 

*■ 

* 

* 

* 

11 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

* 

* 


TABLE V. Damped modes for gravitational perturbations of an extremal Kerr black hole. For a given set of parameters t and 
m, the fundamental frequency is given. Asterisks indicate pairs {I, m) for which no damped modes were found. 


m\£ 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

0 

0.425145 

-0.0718061 

0.665132 

-0.0767351 

0.890510 

-0.0785751 

1.109688 

-0.0794311 

1.325702 

-0.0799001 

1.539879 

-0.0801851 

1.752892 

-0.080372i 

1.965120 

-0.0805011 

2.176794 

-0.0805941 

2.388061 

-0.0806631 

2.599020 

-0.0807161 

1 

0.581433 

-0.0382551 

0.795283 

-0.0589651 

1.013254 

-0.0667191 

1.228731 

-0.0705621 

1.442442 

-0.0728171 

1.655025 

-0.0742871 

1.866861 

-0.0753181 

2.078182 

-0.0760791 

2.289134 

-0.0766631 

2.499811 

-0.0771241 

2.710280 

-0.0774981 

2 

* 

* 

1.028553 

-0.0185721 

1.192475 

-0.0446951 

1.387017 

-0.0560011 

1.589346 

-0.0620971 

1.794694 

-0.0658651 

2.001500 

-0.0684081 

2.209111 

-0.0702351 

2.417208 

-0.0716071 

2.625619 

-0.0726731 

2.834244 

-0.0735261 

3 

X 

* 

1.503222 

-0.0043711 

1.615078 

-0.0313311 

1.783228 

-0.0455081 

1.969595 

-0.0536511 

2.164245 

-0.0588581 

2.363378 
-0.0624501 

2.565215 

-0.0650651 

2.768810 

-0.0670491 

2.973612 

-0.0686041 

4 

X 

X 

* 

* 

* 

2.059795 

-0.0196671 

2.199039 

-0.0356581 

2.367203 

-0.045480i 

2.549284 

-0.0519681 

2.739199 

-0.0565311 

2.933980 

-0.0598971 

3.131996 

-0.0624731 

5 

X 

X 

X 

* 

* 

* 

* 

2.523730 

-0.0101121 

2.632115 

-0.0266681 

2.780307 

-0.0377321 

2.948312 

-0.0453021 

3.127738 

-0.0507321 

3.314388 

-0.0547911 

6 

X 

X 

X 

X 

* 

* 

* 

* 

3.004516 
-0.0030971 

3.080548 

-0.0186731 

3.207368 

-0.0305001 

3.360070 

-0.0389291 

3.527952 

-0.0451101 

7 

X 

X 

X 

X 

X 

* 

* 

* 

* 

* 

3.542772 

-0.0117841 

3.647099 

-0.0238461 

3.783490 

-0.0328981 

8 

X 

X 

X 

X 

X 

X 

* 

* 

* 

* 

* 

4.017480 

0.006130 

4.098411 

-0.0178211 

9 

X 

X 

X 

X 

X 

X 

X 

* 

* 

* 

* 

* 

* 

4.503570 

-0.0019371 

10 

X 

X 

X 

X 

X 

X 

X 

X 

* 

* 

* 

* 

* 

* 

11 

X 

X 

X 

X 

X 

X 

X 

X 

X 

* 

* 

* 

* 

12 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

* 

* 


the parameters s, £, and n, demonstrate that the quasi¬ 
normal frequencies for extremal RN black holes are in 
excellent agreement with the QN frequencies of a near 
extremal black hole for qQ < 0.1. In particular, for s = 0, 
as qQ 0, our results approach the values calculated in 
Ref. [M] for neutral perturbations. As \qQ\ is increased, 
however, we find unexpected results. 

The first unexpected result concerns negative qQ (elec¬ 
tromagnetic attraction). Specifically, as qQ becomes 
more negative, it has been recently found [23l|35] that, for 
non-extremal black holes, there is a special point at which 
the real part of the fundamental quasinormal frequency 


becomes zero and its quasinormal branch is suddenly in¬ 
terrupted. In our numerical simulations for an extremal 
black hole, on the other hand, we have observed that the 
QNM branch does not disappear. This behaviour seems 
to be very general, occurring for the fundamental QNM 
and the first overtone of scalar and Dirac perturbations, 
as shown in Fig. A possible explanation for this dif¬ 
ference is the fact that the analysis of convergence of the 
power series expansion in the continued fraction method 
breaks down when w is a pure negative imaginary num¬ 
ber, which is exactly what happens for near extremal 
black holes at the critical point. In principle, this is 
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TABLE VI. Comparison of the QNM frequencies Mto and the associated separation constants A between a near extremal 
{a/M = 0.999) and an extremal Kerr black hole for m < 0. For each set of parameters s, i, and m < 0, we have calculated the 
least damped mode (n = 0) and the the first overtone (n = 1). 


s 

£ 

m 

n 

a/M = 0.999 

a/M = 1 

0 

1 

-1 

0 

Re(Mtj) lm{Mu) Re(A) Im(A) 

0.239462 -0.093828 1.990310 0.008989 

Re{Muj) lm{Muj) Re(A) Im(A) 

0.239424 -0.093821 1.990294 0.009005 

0 

1 

-1 

1 

0.201620 -0.300295 2.009941 0.024115 

0.201573 -0.300280 2.009963 0.024156 

1 

_! 

1 

I 

1 

0 

0.159112 -0.093385 1.046242 -0.019324 

0.159089 -0.093376 1.046236 -0.019324 

1 

0.116224 -0.311875 1.071028 -0.073697 

0.116178 -0.311858 1.071013 -0.073704 

-1 

1 

-1 

0 

0.204380 -0.091352 2.186087 -0.071537 

0.204349 -0.091348 2.186231 -0.071589 

-1 

1 

-1 

1 

0.158716 -0.299037 2.190971 -0.246458 

0.158677 -0.299029 2.191157 -0.246657 

-2 

2 

-1 

0 

0.343864 -0.083401 4.395467 -0.079325 

0.343862 -0.083384 4.395797 -0.079358 

-2 

2 

-1 

1 

0.316309 -0.257080 4.400545 -0.251561 

0.316305 -0.257024 4.400904 -0.251671 

-2 

2 

-2 

0 

0.291609 -0.088028 4.721185 -0.198122 

0.291553 -0.088026 4.721727 -0.198286 

-2 

2 

-2 

1 

0.250213 -0.276742 4.674323 -0.637396 

0.250146 -0.276738 4.674850 -0.637954 


TABLE VII. The QNM frequencies Mui of an extremal RN black hole for positive values of the electromagnetic interaction 
parameter qQ. 


S 

e 

n 


qQ = 

0.001 



qQ = 

0.01 



qQ = 

^ 0.1 





Q/M 

= 0.999 

Q/M = 1 

Q/M 

= 0.999 

Q/M = 1 

Q/M 

= 0.999 

Q/M = 1 




Re{Muj) 

Im(Mcj) 

Re(McD) 

Im(Afa;) 

Re{Muj) 

Im(Ma;) 

Re(McD) 


Re(Mcj) 


Ke{Muj) 


0 

0 

0 

0.133959 

-0.095843 

0.133959 

-0.095844 

0.138479 

-0.095815 

0.138478 

-0.095816 

0.185412 

-0.093009 

0.185411 

-0.093018 

0 

0 

1 

0.093464 

-0.330652 

0.093465 

-0.330652 

0.097959 

-0.330643 

0.097960 

-0.330644 

0.142213 

-0.329774 

0.142459 

-0.329781 

0 

0 

2 

0.075582 

-0.588328 

0.075581 

-0.588326 

0.080081 

-0.588326 

0.080079 

-0.588323 

0.126173 

-0.589479 

0.124830 

-0.588012 

0 

1 

0 

0.377854 

-0.089529 

0.378142 

-0.089384 

0.382363 

-0.089528 

0.382657 

-0.089379 

0.428778 

-0.089044 

0.429135 

-0.088838 

0 

1 

1 

0.348615 

-0.276440 

0.348680 

-0.276139 

0.353125 

-0.276434 

0.353187 

-0.276128 

0.398835 

-0.275353 

0.398861 

-0.274988 

0 

1 

2 

0.299098 

-0.486659 

0.298958 

-0.486435 

0.303602 

-0.486649 

0.303458 

-0.486425 

0.348584 

-0.485725 

0.348396 

-0.485507 

0 

2 

0 

0.626504 

-0.088912 

0.627073 

-0.088748 

0.631006 

-0.088913 

0.631582 

-0.088746 

0.676875 

-0.088747 

0.677534 

-0.088542 

0 

2 

1 

0.608278 

-0.269526 

0.608671 

-0.269093 

0.612783 

-0.269528 

0.613179 

-0.269087 

0.658498 

-0.269068 

0.658918 

-0.268547 

0 

2 

2 

0.573255 

-0.458754 

0.573373 

-0.458202 

0.577762 

-0.458754 

0.577877 

-0.458195 

0.623207 

-0.458121 

0.623285 

-0.457499 

1 

? 

? 

? 

? 

? 

? 

? 

? 

1 

? 

f 

0 

0.238550 

-0.087811 

0.238682 

-0.087685 

0.243068 

-0.087805 

0.243204 

-0.087672 

0.290216 

-0.086637 

0.290379 

-0.086432 

1 

0.196977 

-0.281368 

0.196896 

-0.281211 

0.201484 

-0.281352 

0.201397 

-0.281194 

0.246685 

-0.279700 

0.246546 

-0.279545 

2 

0.146203 

-0.515626 

0.146077 

-0.515585 

0.150701 

-0.515617 

0.150574 

-0.515578 

0.195347 

-0.514889 

0.195209 

-0.514871 

I 

0 

0.494176 

-0.088402 

0.494613 

-0.088240 

0.498681 

-0.088403 

0.499125 

-0.088237 

0.544802 

-0.088128 

0.545324 

-0.087913 

1 

0.471418 

-0.269611 

0.471647 

-0.269213 

0.475926 

-0.269611 

0.476154 

-0.269205 

0.521728 

-0.268898 

0.521952 

-0.268408 

I 

§ 

§ 

9 

2 

0.429052 

-0.464470 

0.429001 

-0.464039 

0.433560 

-0.464465 

0.433504 

-0.464030 

0.478881 

-0.463620 

0.478768 

-0.463151 

0 

0.745878 

-0.088496 

0.746584 

-0.088327 

0.750378 

-0.088498 

0.751092 

-0.088325 

0.796105 

-0.088385 

0.796906 

-0.088180 

1 

0.730535 

0.267422 

0.731087 

-0.266956 

0.735038 

-0.267425 

0.735594 

-0.266952 

0.780680 

-0.267105 

0.781278 

-0.266553 

2 

0.700602 

-0.452296 

0.700890 

-0.451648 

0.705107 

-0.452298 

0.705394 

-0.451643 

0.750590 

-0.451822 

0.750861 

-0.451087 


also a problem for extremal black holes. However, as we 
increase the number of terms in the continued fraction 
truncations, the real part of the frequency of the QNMs 
seems to approach zero assymptotically as qQ —oo. 


The second unexpected result occurs for positive qQ 
(electromagnetic repulsion). As qQ increases, the quasi¬ 
normal frequencies of a near extremal and an extremal 
black hole behave quite differently. This can be observed 
in Fig.j^for the fundamental frequency and the first over¬ 
tone of scalar and Dirac perturbations. One might argue 
that, since Re{u}) —qQ as qQ increases [331 [33 [SS] , Jo 
tends to become a positive real number and, therefore, 
convergence of the continued fraction method becomes 
compromised (as discussed in Sec. III). While it is true 
that convergence becomes slower as qQ increases, by in¬ 
creasing the number of terms in the continued fractions, 
we were able to obtain accurate results for the range of 
parameters considered in this study. 


C. Mode stability of extremal black holes 


The final section is devoted to an analysis of the sta¬ 
bility of extremal black holes. As explained in Ref. [37], 
there exists three different notions of black hole stability: 
mode stability, linear stability and non-linear stability. 
The first one does not necessarily imply the second one, 
and the second one does not imply the third. Mode sta¬ 
bility is related to the existence of QNMs with positive 
imaginary parts and, therefore, can be studied with the 
method described in this paper. 


Mode stability was established in Ref. [38] for non¬ 
extremal Kerr black holes and in Ref. [33] for non¬ 
extremal Kerr-Newman black holes. In this work we ad¬ 
dress whether or not extremal Kerr and RN black holes 
are mode stable by searching for QNMs with positive 
imaginary parts. For extremal Kerr black holes, we start 
by solving, for several values of s, i and m, the continued 
fraction equations ( 3.11| ) and (3.19) (with 1000 terms) 
using as an initial guess for Moj different values with 
positive imaginary part [for A we use (£ — s)(£-|- s -I-1) as 
the initial guess]. If a solution Mw with positive imagi- 
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TABLE VIII. The QNM frequencies Mui of an extremal RN black hole for negative values of the electromagnetic interaction 
parameter qQ. 


s 

i 

n 


qQ = 

-0.001 



qQ = 

-0.01 



qQ = 

-0.1 





Q/AT 

= 0.999 

Q/M = 1 

07^ 

= 0.999 

Q/M = 1 

Q/AT 

= 0.999 

Q/M = 1 




^e(Muj) 

Im(Mcj) 

Re(Mci;) 

Im(A7a;) 

Re{Muj) 

Im(Ma;) 


Im(Mu;) 

Re(Mcj) 


Re(Mu;) 

lm{Mio) 

0 

0 

0 

0.132959 

-0.095843 

0.132959 

-0.095844 

0.128479 

-0.095815 

0.128478 

-0.095816 

0.085413 

-0.093016 

0.085411 

-0.093018 

0 

0 

1 

0.092465 

-0.330652 

0.092465 

-0.330652 

0.087960 

-0.330643 

0.087960 

-0.330644 

0.042459 

-0.329779 

0.042459 

-0.329781 

0 

0 

2 

0.074582 

-0.588327 

0.074581 

-0.588326 

0.070079 

-0.588324 

0.070079 

-0.588323 

0.024829 

-0.588012 

0.024830 

-0.588012 

0 

1 

0 

0.376855 

-0.089528 

0.377142 

-0.089384 

0.372376 

-0.089518 

0.372657 

-0.089379 

0.328906 

-0.088937 

0.329135 

-0.088838 

0 

1 

1 

0.347615 

-0.276439 

0.347680 

-0.276139 

0.343119 

-0.276422 

0.343187 

-0.276128 

0.298773 

-0.275229 

0.298861 

-0.274988 

0 

1 

2 

0.298097 

-0.486659 

0.297958 

-0.486435 

0.293592 

-0.486649 

0.293458 

-0.486425 

0.248492 

-0.485727 

0.248396 

-0.485507 

0 

2 

0 

0.625506 

-0.088911 

0.626073 

-0.088748 

0.621023 

-0.088906 

0.621582 

-0.088746 

0.577046 

-0.088671 

0.577534 

-0.088542 

0 

2 

1 

0.607278 

-0.269524 

0.607671 

-0.269093 

0.602788 

-0.269511 

0.603179 

-0.269087 

0.558553 

-0.268902 

0.558918 

-0.268547 

0 

2 

2 

0.572254 

-0.458753 

0.572373 

-0.458202 

0.567755 

-0.458740 

0.567877 

-0.458195 

0.523138 

-0.457984 

0.523285 

-0.457499 

1 

? 

? 

? 

? 

? 

? 

? 

? 

2 

1 

? 

? 

0 

0.237550 

-0.087810 

0.237682 

-0.087685 

0.233075 

-0.087792 

0.233204 

-0.087672 

0.190278 

-0.086505 

0.190379 

-0.086432 

1 

0.195976 

-0.281368 

0.195896 

-0.281211 

0.191473 

-0.281350 

0.191397 

-0.281194 

0.146583 

-0.279689 

0.146546 

-0.279545 

2 

0.145203 

-0.515627 

0.145077 

-0.515585 

0.140698 

-0.515621 

0.140574 

-0.515578 

0.095320 

-0.514930 

0.095209 

-0.514871 

I 

0 

0.493178 

-0.088401 

0.493613 

-0.088240 

0.488697 

-0.088394 

0.489125 

-0.088237 

0.444962 

-0.088034 

0.445324 

-0.087913 

1 

0.470418 

-0.269610 

0.470647 

-0.269213 

0.465926 

-0.269594 

0.466154 

-0.269205 

0.421728 

-0.268728 

0.421952 

-0.268408 

I 

§ 

§ 

2 

0.428051 

-0.464469 

0.428001 

-0.464039 

0.423549 

-0.464457 

0.423504 

-0.464030 

0.378770 

-0.463539 

0.378768 

-0.463151 

0 

0.744880 

-0.088495 

0.745584 

-0.088327 

0.740396 

-0.088491 

0.741092 

-0.088325 

0.696286 

-0.088319 

0.696906 

-0.088180 

1 

0.729536 

-0.267420 

0.730087 

-0.266956 

0.725047 

-0.267409 

0.725594 

-0.266952 

0.680773 

-0.266944 

0.681278 

-0.266553 

2 

0.699601 

-0.452294 

0.699890 

-0.451648 

0.695105 

-0.452281 

0.695394 

-0.451643 

0.650564 

-0.451653 

0.650861 

-0.451087 


10 

9 

8 

7 

6 

E 5 
4 
3 
2 
1 
0 


FIG. 4. (Colors online.) Existence of DMs is summarized in 
this plot. Pairs (i, m) for which no DMs are allowed are rep¬ 
resented by a black cross. Red triangles indicate the existence 
of DMs for s = —1/2. Blue dots indicate the existence of DMs 
for all bosons (s = 0, 1 and 2 simultaneously), magenta dots 
indicate the existence of DMs for s = 0 and 1, but no DMs 
for s = 2, and, finally, orange dots indicate the existence of 
DMs only for s = 0. The dashed line corresponds to the the¬ 
oretical result m/(^-|-1/2) = 0.744 obtained using the eikonal 
approximation. 



tremal Kerr black holes nor for extremal RN black holes), 
indicating that these extremal black holes, like their non¬ 
extremal counterparts, are mode stable. 

In order to illustrate our search, we exhibit some of 
the plots obtained for scalar perturbations with £ = 1 
and qQ = 10 around an extremal RN black hole. In 
particular, note that the plots in Fig. suggest that so¬ 
lutions with positive imaginary parts might exist around 
Moj = 0 and around Mlo = qQ = 10. However, as the 
number of terms in the continued fractions is increased, 
the minimum points of the contour plots move towards 
Im(Ma;) = 0 (see Fig. [^. The same behaviour is found 
when attempting to solve equation (3.11): as the num¬ 
ber of terms in the continued fractions is increased, its 
solutions move towards the real axis, meaning that the 
continued fractions do not converge. 

Conversely, the (stable) quasinormal modes found us¬ 
ing our method are indeed convergent solutions. In Fig.[^ 
we illustrate this result graphically for the case s = 0, 
i = 1, qQ = 0.1. One can see that, as the number of 
terms in the continued fractions is increased, the loca¬ 
tion of the maxima remain unchanged. As explained in 
Sec. |IV[ every time a solution was foun d, we increased 
the number of terms in equation (3.11) to confirm the 
result. 


nary part is found, the procedure is repeated with more 
terms in the continued fractions in order to confirm the 
results. For extremal RN black holes we proceed simi¬ 
larly, solving the continued fraction equation (3.11) for 


several values of s, i, and qQ. Additionally, we also make 
contour plots of the logarithm of the LHS of equation 
(3.11) as a function of Re(Ma;) and Im(Ma;). If any 


unstable mode is suspected from such plots, we use the 
corresponding value of the frequency as an initial guess 
in the root finding algorithm to confirm the result. 

We report that no quasinormal frequencies with pos¬ 
itive imaginary parts have been found (neither for ex- 


V. FINAL REMARKS 

In this paper, inspired by the ideas of Ref. m, we have 
successfully implemented a continued fraction method to 
determine the QNMs of neutral massless perturbations 
around an extremal Kerr black hole and the QNMs of 
charged massless perturbations around an extremal RN 
black hole (we remark the unified framework employed, 
which is applicable in both situatio ns). Starting with the 
perturbation equations (2.2) and (2.3), we set (exactly) 


a = M OY Q = M and, using the asymptotic behaviour 
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FIG. 5. (Colors online.) Real and imaginary parts of the QNMs for negative values of qQ. The blue (dotted) and the red 
(dashed) curves correspond, respectively, to scalar {s = j = 0) and Dirac (s = —1/2, i = 1/2) fields around a near extremal 
(QIM = 0.999) RN black hole, while the magenta (solid) and the black (dot-dashed) curves correspond to scalar (s = £ = 0) 
and Dirac (s = —1/2, £ = 1/2) fields around an extremal RN black hole. The top and bottom plots correspond, respectively, 
to fundamental QNMs and its first overtones. Note the critical value of qQ below which the quasinormal branches disappear 
for near extremal black holes. 


given in (12.51) and (2.6), write a power series solution of 


the radial equation around the ordinary point r = 2M. 

We were able to obtain the quasinormal frequencies 
for arbitrary values of the parameters s, £, m and qQ. 
Our results agree with Leaver’s original continued frac¬ 
tion method for near extremal black holes, as explicitly 
shown in tables [Ilvllvinl Our results also agree with 
previous results in the literature for near extremal Kerr 
black holes miMiiin] and for near extremal RN black 
holes |23l 135). In particular, for m ^ 0, our numerical 
calculations, as shown in tables |n| and [Vl and in Fig. 
are compatible with the scalar and gravitational damped 
modes of near extremal Kerr black holes, which were 
studied in Refs. [301 [SI]. Besides extending the analy¬ 
sis of Refs. [Sni ISI] to Dirac and electromagnetic fields, 
we have studied in detail the behaviour of the most stable 
damped modes as the parameters I and m are varied. 

It is important to note that the asymptotic behaviour 
of the solution of the radial perturbation equation of a 
near extremal black hole is not the same as (2.5), which 


holds only for extremal black holes. The remarkable 


agreement that we have observed between the QNMs of 
near extremal and extremal black holes, for most sets of 
parameters s,£, m and gQ, supports the idea that the 
QNM limit of extremal black holes is continuous. Never¬ 
theless, there are some important differences, as observed 
in Figs. [^ and [^ In particular, contrary to extremal RN 
black holes, near extremal black holes exhibit critical val¬ 
ues of the electromagnetic parameter below which the 
quasinormal branches disappear. Moreover, as the elec¬ 
tromagnetic interaction increases, notable differences be¬ 
tween the quasinormal frequencies of extremal and near 
extremal RN black holes emerge. 

We would also like to remark that another commonly 
used method, the WKB method, besides being an ap¬ 
proximation, becomes less accurate as the extremal limit 
is approached [H] (see also Ref. [31] for an implementa¬ 
tion which works for near extremal black holes, but only 
in the eikonal limit). Since the original continued fraction 
method fails when a = M oy Q — M, our implementation 
of the ideas of Ref. m in such a case provides the most 
accurate method to determine the QNMs of an extremal 



















13 



FIG. 6. (Colors online.) Real and imaginary parts of the quasinormal frequencies for positive value of qQ. The blue (dotted) 
and the red (dashed) curves correspond, respectively, to scalar [s = j = 0) and Dirac (s = —1/2, i = 1/2) fields around a 
near extremal (QjM = 0.999) RN black hole, while the magenta (solid) and the black (dot-dashed) curves correspond to scalar 
(s = ^ = 0) and Dirac (s = —1/2, i = 1/2) fields around an extremal RN black hole. The top and bottom plots correspond, 
respectively, to fundamental QNMs and its first overtones. Note the difference between extremal and near extremal black holes 
as the electromagnetic interaction qQ is increased. 


Kerr black hole and the QNMs of charged perturbations 
around an extremal RN black hole. 

Finally, we point out that we have found no quasinor¬ 
mal frequencies which possess positive imaginary parts, 
indicating that extremal Kerr and extremal RN black 
holes are mode stable. Onozawa et. al, similarly, did 
not report finding unstable modes related to (uncharged) 
perturbations around an extremal RN black hole El- 
Recently, however, extremal RN and Kerr black holes 
have been shown to be linearly unstable. The first re¬ 
sults were obtained by Aretakis in Refs. [42ti4^ . Later, 
inspired by Aretakis’ works, Lucietti and Reall have 
shown that extremal Kerr black holes are unstable un¬ 
der linearized gravitational and electromagnetic pertur¬ 
bations. There might be, however a connection between 
zero-damping modes and these linear instabilities. While 
for near extremal black holes the zero-damping modes 
have small (but non-zero) negative imaginary parts, it is 
believed that, for extremal black holes, they are purely 
real [ISHIS] • There has been some suggestions that these 
purely real modes are related to instabilities [531 US] , but 


whether this is really the case (and their relation with 
these linear instabilities) has yet to be established. 
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FIG. 7. (Colors online.) Searching for mode instabilities of 
the extremal RN black hole for s = 0, ^ = 1, and qQ = 10: 
contour plots of the logarithm of the LHS of equation (3.111 
as a function of Re(Ma;) and Im(Ma;). The number of terms 
in the continued fractions is 1000. The bottom plot is a close- 
up view of the upper one and suggests that unstable modes 
with small \m{Muj) might exist. 


Appendix: Recurrence relation, gaussian elimination 
and decoupling 


The coefficients for the 5-term recurrence relation (3.4) 
are given explicitly by 


an=v? + n, Pn = Pin, (A.l) 

7„ =-2n^-I-P2n-I-/a, 5n = PAn + P^, (A.2) 

e„ = + Pen + P 7 , (A.3) 


where 

Pi = 4:{im — iMco — s), (A.4) 

P 2 = 2{1-Aim+mMuj), (A.5) 

P 3 = - 2[1 -H 2s -b 2A -b 2iMuj{A + 23iMui) (A. 6 ) 
— 2zto(1 — 6 iMuj)], (A.7) 

P 4 = A{im — iMuj + s), (A. 8 ) 

Pe =—4{im — iMuj + s){l + AiMuj), (A.9) 

Pg = -3 - 8 iMu}, (A.IO) 

P 7 = 2[l-b6iMa;-8(Mw)^], (A.ll) 

for an extremal Kerr black hole, and 

Pi=AiiqQ-s), (A.12) 

P 2 = 2{l-6iqQ + l2iMuj), (A.13) 

P 3 =-2[l-31gQ-b2s-b2A-b6*Ma; (A.14) 

- 2{qQ - AMu;){3qQ - AMu})], (A.15) 

P 4 = A{iqQ + s), (A.16) 

Pg = -A{iqQ + s)(l - 2iqQ + AiMu), (A.17) 

Pe = -3 + AiqQ-8iMuj, (A.18) 

P 7 = 2(1 - 2iqQ + AiMuj){l - iqQ + 2iMuj), (A.19) 


for an extremal RN black hole. In order to reduce the 


5-term recurrence relation (3.4) to a 3-term recurrence 


relation, we use a double gaussian elimination proce¬ 
dure BM- We define 


ai = q;i, P[ = Pi, = 71, 

0^2 = 012, P2 = P2, 72 = 72, ^2 = <^2 

= 0, aj, = an, P'^ = Pn - 

In — In — JP~P'n-l,5'n = — gi " 


followed by 


n > 3 
n > 3 
(A.20) 


a" = a'l, P'' = P[, = 7'i, 

P' = a' P" = P' _^ 

u / 

7n = 7n - 




-a 


n—l'> 


'' — n 


n>2 (A.21) 

n> 2 


transforming eqs. (3.2)-(3.4| into 


o'n^n+l + Pn^n + Tn^n-l — 0, 71 > 1. (A.22) 

This 3-term recurrence relation couples even {a 2 n) ^^nd 
odd (a 2 n-i-i) expansion coefficients. In order to deco uple 
them, we proceed as follows. First, we write (A.22) for 


277-1-1 and solve for a 2 r!,-i- 2 - Next, we write (A.22) for 


2 n, solve for a 2 n, and plug the result into the expression 
obtained earlier for a 2 n-i- 2 , resulting in 

'y2n+li^2nO'2n+l + l2n^2n-l) 


<l2n+2 — 


"'2n-H 

'^'in+l 


n 2 n+l 


^2n+l/^2n 


P'^ 


(A.23) 


Finally, we substitute this expression into ( A.22| ) for 
2 n + 2 , successfully decoupling the odd terms from the 


even: 


On^n+l A~ Pn^n A~ Pn^n—l — 0 , 


77 > 1, (A.24) 
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FIG. 8. (Colors online.) Searching for mode instabilities of the extremal RN black hole for s = 0, £ = 1, and qQ = 10: close-up 
view of the contonr plots in Fig. [^aronnd Muj = 0 (top panels) and around Mlo = qQ = 10 (bottom panels). From left to 
right, the panels correspond respectively to 1000, 1400 and 2000 terms in the continued fractions. As the number of terms is 
increased the minima move towards Im(Ma;) = 0, where we cannot guarantee converge of the continued fractions (as explained 
in Sec. |III[|. This indicates that no unstable mode exists around Mlo = 0 or Moo = 10. 



FIG. 9. (Colors online.) Contour plots indicating the location of the first three damped modes of the extremal RN black hole 
for s = 0, £ = 1, and qQ = 0.1. From left to right, the panels correspond respectively to 1000, 1400 and 2000 terms in the 
continued fractions. Unlike the contour plots in Fig. as the number of terms is increased, the minima do not move. The 
precise location of these quasinormal modes is calculated in Table [VTIl We also note that the point located at Mlo ~ 0.38 —0.3i 
is a local maxima and, consequently, not a solution of < |01 ). 


where dn = cL 2 n+i and the recurrence coefficients are Analogously, it is possible to write a 3-term recurrence 
given by relation for the even terms: 


a 


O 

n 


P 


o 

n 


— ^2n+2t 

_ „// /^2n + 1^2n.+2 

— l2n+2 

a" " n 

— n" B" 

“ 271 + 1^271 



(A.25) 


= 0 , 


In 


Pn^n 


n > 1, (A.26) 
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where c„ = a 2 n and 


— 0 ^ 271+11 

Pe _ ,,// l^2n^2n+l , ^2n + 1^2n - l'y2r 

Pn - 72n+l a''^ + «L/32„-l 

e _ /^ 2 ji + l72ra-l72Ti 

in ~ , 


(A.27) 


Note that, while the original decoupling procedure of 
Ref. [14] only works when i5„ = 0, our decoupling works 
for Sn 0 (which is the case in all scenarios we have 
analysed in this paper). 
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